Spin waves in the Frustrated Kagome Lattice Antiferromagnet KFe 3 (OH) 6 (S04)2 
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The spin wave excitations of the ideal S = 5/2 Kagome lattice antiferromagnet KFe3(OH)e(S04)2 
have been measured using high-resolution inelastic neutron scattering. We directly observe a flat 
mode which corresponds to a lifted "zero energy mode," verifying a fundamental prediction for the 
Kagome lattice. A simple Heisenberg spin Hamiltonian provides an excellent fit to our spin wave 
data. The antisymmetric Dzyaloshinskii-Moriya interaction is the primary source of anisotropy and 
explains the low temperature magnetization and spin structure. 
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Geometrically frustrated spin systems have received 
considerable attention in recent years due to the pres- 
ence of remarkable properties such as spin iceQJ^], spin 
nematicQ, and spin liquid behaviors H uW The 
Kagome lattice antiferromagnet is a highly frustrated 
two-dimensional lattice, being comprised of corner- 
sharing triangles. For classical Heisenberg spins, the 
ground state of a Kagome antiferromagnet is infinitely 
degenerate, but the spins are believed to order in the 
T — > limit by a process known as "order by disorder" 
7, 8]. On the other hand, there are predictions that the 
ground state of the S = 1/2 Kagome lattice is disordered, 
being a realization of the long sought after quantum 
spin-liquid 0, 0, 0, 0]. Experimentally, several materials 
have been studied that are believed to be realizations of 
the Kagome lattice antiferro magnet , such as scgo|I3, 
volborthitefll). and iarosites |l2l llH, . However, these 
materials are often plagued by non-stoichiometry issues 
or have structural differences from the ideal Kagome net- 
work. In this paper, we present a high-resolution neutron 
scattering study on a single crystal sample of a pure sys- 
tem with an ideal Kagome net work |l5j. the iron jarosite 
KFe 3 (OH) 6 (S04)2- This allows us to directly compare 
our data with fundamental theoretical predictions. 

One of the hallmarks of highly frustrated systems is 
the presence of "zero energy modes" which result from 
the huge ground-state degeneracy. For the Kagome lat- 
tice Heisenberg model, the only constraint for the ground 
state is that the spins on each triangle be oriented 120° 
relative to each other. A "zero energy mode" for the 
Kagome lattice is depicted in Fig. 1(a). The small loops 
at the tips of the spins illustrate rotations of two of the 
spin sublattices about the axis defined by the third spin 
sublattice. These spins, forming a chain, can collectively 
rotate around the loop paths with no change in energy 
(the 120° angles on each triangle are maintained). Fur- 
thermore, the spins on different parallel chains can be 
excited independently. Hence, this t ype of excitation 
costs no energy and is non-dispersive [lfjj. This mode 



has not been directly observed previously, and, since 
it occurs at zero energy, it is difficult to measure with 
most experimental techniques. Here, we report a di- 
rect observation of a zero energy mode in iron jarosite, 
which is lifted to a finite energy due to the presence of 
spin anisotropy resul ting from the Dzyaloshinskii-Moriya 
(DM) interaction [l7l[l8| . 

KFe3(OH)g(S04)2 is an ideal realization of a Kagome 
Heisenberg antiferromagnet due to its undistorted lat- 
tice, full y occup ied magnetic sites, and weak interlayer 
coupling |l5j, LUj. The 5 = 5/2 Fe 3+ ions form a perfect 
Kagome arrangement and are surrounded by an octahe- 
dral environment of oxyg ens. Previous powder neutron 
diffraction studies [T3. Il3l| identified the long-range mag- 
netic order below Tjy = 65 K as that shown in Fig. 1(a). 
For this spin ar rang ement, each triangle has positive vec- 
tor chirality^, llpJ ]. such that the spins rotate clockwise 
as one traverses the vertices of a triangle clockwise. Mag- 
netization results on single crystals have determined that 
the spins are canted slightly out of the Kagome plane, 
yielding an "umbrella" structure. Moreover, each plane 
has a net ferromagnetic moment, with an antiferromag- 
netic coupling between adjacent planes [l9|. 

We first studied the magnetic excitations using a 
deuterated powder sample (mass = 4.92 g) on the DCS 
and BT4 spectrometers at the NIST Center for Neutron 
Research, as shown in Figures 1(b) and 1(c). Figure 1(c) 
shows a difference plot of the intensity as a function of 
neutron energy loss measured above (T = 70 K) and be- 
low (T = 13 K) the Neel temperature. This difference 
plot removes most of the phonon contributions to the 
spectrum, hence yielding the spin wave density of states. 
Despite the powder average, the spectrum shows one 
sharp feature at hu>o ~ 8 meV and a second broad peak at 
about 2luq- Both features appear as excitation bands over 
a wide range of \Q\, as shown in Fig. 1(b). This behavior 
is quite similar to that observed in strongly frustrated 
spinel systems where the excitation at u>o has been de- 
scribed as a local resonance [ioj. It is tempting to identify 
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FIG. 1: (color online) (a) The ground state spin configuration 
of iron jarosite. The magnetic unit cell, shown by the yellow 
shaded area, is the same as the chemical unit cell. The dotted 
loops illustrate the zero energy excitations as described in the 
text. The Cartesian axes are used in the discussion of the DM 
interaction, (b) Intensity contour map of the inelastic scat- 
tering spectrum at T = 4 K of a powder sample measured 
using the time-of-flight DCS spectrometer with an incident 
neutron wavelength of 1.8 A. (c) Inelastic neutron scattering 
measured on a powder using the BT4 spectrometer with col- 
limations 40' — 20'. The data show the difference between the 
intensities above and below the Neel temperature Tjv = 65 K 
and is a measure of the spin wave density of states. 



the features observed in the excitation spectrum as one- 
and two-magnon scattering since strong multi-magnon 
scattering might be expected due to the strong frustra- 
tion and cubic terms in the spin-Hamiltonian. However, 
as shown below, our single crystal measurements pro- 
vide much greater detail and demonstrate that these are 
regular spin wave modes, albeit with unusual dispersive 
behavior. 

The spin wave dispersions were obtained from inelas- 
tic neutron scattering measurements on a single crys- 
tal sample (composed of four co-aligned crystals of to- 
tal mass 101 mg) grown using a hydrothcrmal method 
reported previously 0, 0|. High-resolution measure- 
ments were performed using the triple-axis spectrom- 
eter HB1 at the High Flux Isotope Reactor at Oak 
Ridge National Laboratory with the sample aligned in 
the (HKO) and (HHL) zones with the final energy fixed 
at either 13.6 meV or 14.7 meV. Vertically focused py- 
rolytic graphite (PG) crystals were used to monochro- 
mate and analyze the incident and scattered neutron 
beams using the (0 2) reflection. Horizontal collima- 
tions of 48' — 60' — sample— 40' — 120' were employed, and 
PG filters were placed in the scattered beam to reduce 
higher-order contamination. The sample was cooled to 
T = 10 K using a closed cycle 4 He cryostat. 
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FIG. 2: (color online) (a) Energy scans at Q = (1 0) and 
(1.1 0). (b) Q-scan at hui = 9.5 meV. (c) Energy scans at 
Q = (1 1 0), (1.25 1 0), and (1.5 1 0). The solid lines show the 
fits to the spin wave dispersion relation described in the text, 
convoluted with the instrumental resolution function. In the 
lower panel of (c), the CF prediction is shown by the red line 
and the DM prediction by the blue line. 

A series of energy scans (at constant Q) and Q-scans 
(at constant energy) were performed, and a few represen- 
tative scans are shown in Fig. 2. The observed peaks were 
initially fit with narrow Gaussians convoluted with the 
experimental resolution function. Subsequent fits were 
performed taking into account the empirical dispersion 
of the excitations. The peaks are resolution-limited, and 
the line-shapes are simply governed by the convolution 
with the instrumental resolution. A summary of all of 
the peak positions and intensities is shown in Figure 3(a). 
The error bars plotted in Figs. 3(a) and (b) correspond 
to three times the statistical error or one-tenth of the 
instrumental resolution, whichever is larger. The most 
striking feature of the data is the relatively flat mode 
near 7 meV which barely disperses, even out to the zone 
boundary. 

The energy scan in Fig. 2(a) at the magnetic Brillouin 
zone center (1 0) reveals two spin gaps, one at 1.8(1) 
meV (which is non-degenerate), and the other at 6.7(1) 
meV (which is two-fold degenerate within the experimen- 
tal resolution). At Q = (1.1 0), the lower-energy mode 
has dispersed to higher energy and merges with the flat 
mode located around 7 meV. The other upper-energy 
mode disperses strongly, moving to a high zone boundary 
energy of about 19 meV. Figure 2(c) shows constant-Q 
scans of the flat mode within a Brillouin zone centered at 
(1 1 0). This excitation barely disperses, starting from 
about 7 meV at the zone center and reaching about 9 
meV at the zone boundary. We identify this flat mode 
as the "zero energy mode" of the Kagome lattice which 
is lifted in energy for reasons discussed below. 

We have fit the observed spin wave dispersions using 
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FIG. 3: (color online) Spin wave dispersion along the high 
symmetry directions in the 2D Brillouin zone. As discussed 
in the text, the solid lines denote a fit to the DM model (in a) 
and to the CF model (in b). The model parameters are given 
in units of meV. (c) Wave vector dependence of the spin wave 
intensities. The solid lines correspond to (n[ui) + l)/w(Q), 
where n(u>) is the Bose occupation factor and ui(Q) is obtained 
from the DM model. The data take into account the de- 
convolution with the instrumental resolution and the Fe 3+ 
magnetic form factor. 



the following generic Hamiltonian: [2 lj 
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where S nTl (£ nml ) indicates summation over pairs of 
nearest neighbors (next nearest neighbors), — 
[0, D y (i,j), D z (i,j)] is the DM vector for bond i — j as 
shown in Fig. 1(a), and the single- ion anisotropy terms 
are those used by Nishiyama et al. Ti\ in their treatment 
of the spin wave spectrum in jarosites. Here the primed 
spin components refer to the local axis associated with 
the rotated oxygen octahedra (see Ref. for details). 
We ignore the weak interplane coupling, which is several 
hundred times smaller than J\ [19| ■ 

The DM interaction is allowed for this crystal struc- 
ture and merits further discussion J^. For bond 1- 
2, it has the form D 1]2 = (0,D y ,D z ). The other 
DM vectors can be obtained from symmetry, such as 
D,(l,2) = D*(2,3) = D,(3,l), D„(4,5) = -D„(l,2) 
and £> 2 (4, 5) = — D z (l,2). Note that the direction of 
the DM vector oscillates from bond to bond along the x- 
direction. The z-component of the DM vector favors the 
spins to lie in the ab— plane and therefore effectively acts 



like an easy-plane anisotropy. The sign of D z breaks the 
symmetry between positive and negative vector chirality. 
Similarly the D y component breaks the rotational sym- 
metry around the c-axis and creates a small anisotropy 
with respect to in-plane orientations. The effect of D y is 
also to cant the spins so that they have a small out-of- 
plane component, consistent with the observed umbrella 
spin configuration. 

We may describe the spin wave data in terms of two 
simple spin models. In the first of these, which we call 
the DM model, we neglect the single ion anisotropy, so 
that the only nonzero parameters are J\, J2, D y , and 
D z . In the second model, which we call the CF (crystal 
field) model, all the anisotropy is attributed to the single- 
ion crystal field, so that the only nonzero parameters are 
Ji, J21 D, and E. In both cases, J\ is the dominant 
interaction. The numerical results obtained from these 
two models are plotted as the lines in Fig. 3, and the 
approximate analytic expressions for the spin gaps at the 
r point are given in Table [I] 

To account for the observed umbrella spin structure, 
we considered the effect of spin-canting on the spin wave 
energies. We find that the splitting of the mode energies 
at the high symmetry points arc particularly sensitive 
to the magnitude of the spin canting angle out of the 
Kagome plane. The best fit to the DM model is depicted 
by the lines in Fig. 3(a) and describes the data very well. 
We reproduce not only the gaps at the zone center, but 
also the small dispersion of the flat-mode. This small 
dispersion is a result of a weak nnn interaction J2- Wc 
note that J2 is positive (antiferromagnetic), which favors 
the experimentally observed ground state. The D z com- 
ponent of the DM vector also reinforces selection of this 
ground state. The "zero energy mode" is lifted by an en- 
ergy equal to the out-of-plane spin wave gap, consistent 
with the spin rotations depicted in Fig. 1(a). The gaps 
at the r point obtained numerically are in good agree- 
ment with the analytic results given in Table [I] The 
DM model yields a spin-canting angle of about 2.1°, in 
reasonable agreement with the experimentally estimated 
value of 1.5(5)° 0. 

The lines in Fig. 3(b) show the best fit to the CF model, 
which also is in reasonable agreement with the data. This 
model predicts a canting angle of 0.8°, which is slightly 
smaller than, but consistent with, the deduced experi- 
mental value. However, the least-squares goodness-of-fit 
to the spin wave spectrum is not as good as that us- 
ing the DM model. The difference is most apparent in 
the numerical results for the mode splitting at the zone 
center. The CF model yields a relatively large splitting 
of about 0.71 meV for the 7 meV mode at the T point, 
whereas the data indicate that this splitting is smaller 
than 0.4 meV. The red line in Fig. 2(c) shows the CF 
prediction for the energy scan at the zone center, and we 
see that this does not describe the lineshape very well. 
The DM model, depicted by the blue line, describes the 
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TABLE I: Spin wave energies at the zone center for the DM 
model and the CF model. Here, Co = uj/S, J = Ji + J2, and 
8 ~ 20° is the oxygen octahedra tilting angle. 
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data better. Moreover, as pointed out in Ref. |23j], the 
single ion anisotropy of the Fe 3+ ion is expected to be 
small since it appears at second order in the spin-orbit 
coupling, whereas the DM term appears at first order. 

Therefore, we believe the observed spin wave spectrum 
is most naturally explained by a simple model which has 
only nearest and next-nearest isotropic interactions plus 
the DM interaction. The obtained fit parameters (in 
meV) are Ji = 3.23(4), J 2 = 0.11(1), \D y \ = 0.218(5), 
and D z = —0.195(5). From a previous study[l9j, a value 
for J\ of 3.9(2) meV was obtained from a fit of the suscep- 
tibility to a high-temperature series expansion resultjlfij. 
The values of J\ are in reasonable agreement, and the 
agreement would be even closer if the effects of J2 and 
the DM term were taken into account in the suscepti- 
bility fit. As a further comparison, the susceptibility 
in Ref. indicated a value of Ag/g ~ 0.06, where 
g is the free electron Lande factor and Ag is its shift 
in the crystalline environment. The magnitude of the 
DM vector can be estimated from Moriya's calculation 
as \D\/Ji ~ Ag/g From the current study, we have 
Dy/Ji ~ \D Z \/J\ ~ 0.06 ~ Ag/g, showing very good 
agreement between measurements of the spin dynamics 
and the bulk thermodynamics. 

Finally, from the analytic expressions for the spin gaps 
given in Table [I] we note that the in-plane gap is pro- 
portional to \D y \ while the out-of-plane gaps are propor- 
tional to \J J\ D z . Since J\ is large compared to other 
interactions, the out-of-plane gap is significantly higher 
than the in-plane gap, despite the similar magnitude of 
D y and D z . This also suggests that at high temper- 
atures (even above T N ), the spins would feel an easy- 
plane anisotropy and therefore display ATK-hke spin dy- 
namics at high temperatures. This picture is consistent 
with the previous neutron scattering measurements of 
the critical fluctuations above T/v which indeed have XY 
symmetry That previous study also showed that vec- 
tor chiral order is apparent above Tjy, consistent with 
the presence of the DM term and the positive sign of J2. 
Another interesting aspect of the analytic results is that 
the spin wave spectrum does not depend on the sign of 
D y , which determines the direction of the spin-canting 
relative to the in-plane order. This information on the 
canting cannot be determined from the currently avail- 
able data and would be important for testing microscopic 
calculations of the DM interactions. 



In summary, the spin wave spectrum of a Kagome lat- 
tice antiferromagnet has been measured using inelastic 
neutron scattering. We observe a flat, lifted "zero en- 
ergy mode" at ~7 meV, whose presence reflects the huge 
ground-state degeneracy of the ideal Kagome Heisenberg 
antiferromagnet. We have also determined all of the rel- 
evant spin Hamiltonian parameters by fitting our data to 
a Heisenberg model with the antisymmetric DM interac- 
tion. This realization of the Kagome antiferromagnet is 
perhaps the best characterized geometrically frustrated 
spin system, and, as such, would be useful for precise 
tests of theoretical predictions. These results also high- 
light the importance of single crystal measurements to 
accurately interpret the data taken on powder samples 
of frustrated magnets. 
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